Charge Fluctuations in vr+p and K+p Collisions at 250 GeV/c 
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We report on measurements of event-by-event charge fluctuations in 7r"^p and K^p collisions at 
250 GeV/c. The dependence of these fluctuations on the size of the rapidity windows are presented 
for the first time in the full phase space domain. The corrections for the influence of global charge 
conservation and leading-particle stopping are tested by the data. The discrepancy due to incomplete 
correction given by STAR and PHENIX are estimated. The dependence of the fiuctuations on the 
position of the rapidity bin and on the multiplicity at different rapidity windows are also presented. 
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The subject of event-by-event fluctuations has cur- 
rently drawn a lot of attention in both theoretical 
and experimental studies of relativistic heavy ion colli- 
sions [llQl- It is argued that information on the QCD 
phase transition — the formation of Quark Gluon Plasma 
(QGP) — can be inferred from measurements ^] among 
which event-by-event charge fluctuations are considered 
as a promising signature j^l^. Due to the fractional elec- 
tric charges of quarks, the charges spread more evenly 
throughout the QGP volume than in a hadronic gas and, 
therefore, the fluctuations are expected to suffer an ob- 
servable suppression in a QGP jj]. 

Recently, it has been demonstrated that event-by- 
event charge fluctuation can be directly related to a ther- 
modynamic signature — the anomalous proton-number 
fluctuation at the critical point which is supposed to 
enhance the charge fluctuations. The observed enhance- 
ment of charge fluctuations at RHIC and SPS 0] seems 
to be a good support for these arguments, though the ef- 
fects of limited detector acceptance and other corrections 
need to be further investigated. 

The charge fluctuations are also sensitive to other ef- 
fects, as the number of resonances at chemical freeze- 
out 0, Q and fluctuations occurring in the initial 
stage Q. The corresponding analyses are interesting by 
their own beyond the QGP hypothesis 0- 

There are mainly two kinds of measures for the event- 
by-event charge fluctuations on the market at present, 
others being related to these under reasonable assump- 
tions 2]. One is that of net c/iarge fluctuations, the other 
that of charge ratio fluctuations. The direct measure of 
net charge fluctuations is the variance of net charge Q, 



where n+ and are the numbers of positive and nega- 
tive particles observed in a particular phase space window 
under consideration. The average is over all events in the 
sample. If charge is randomly assigned to each particle, 
SQ^ = (rich), where rich = + n~ . So the measure for 
net charge fluctuations is defined as 



D{Q) 



(rich) 
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i.e., equal to 4 for independent emission. 

In order to reduce the fiuctuations of ri-|_ and n_ 
due to the variation of impact parameter, charge ratio 
R = ri^ /n~ fluctuations are recommended in I3| and the 
corresponding measure is 



D{R) = (rich) • 



(3) 



where 5R'^ = (i?^) - (i?)^. 

In the high multiplicity limit, the above two measures 
are approximately equal, with the leading order correc- 
tion being ~ l/(nch)- 

In accounting for the charge conservation in a large 
rapidity window and a non-zero net charge due to non- 
negligible baryon stopping, two correction factors [l0| . 
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are applied to the D-measures of Eq's. (2) and (3): 
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FIG. 1; (a) Dependence of D{Q) on the size of the central rapidity window |y| < Sy/2 with O.OOlGeV/c < pt < lOGeV/c 
(open circles), pt > O.lGeV/c (open squares), pt > 0.2GeV/c (crosses), and pt > 0.2GeV/c, A<j) = n/2 (open triangles), 
(b) Corrected version of (a) by two factors, (c) Corrected version of (a) by only the Cy factor of Eq. (4). (d) Dependence of 
D{Q) on the position of a unit-width rapidity window. Dotted, dashed, solid and dot-dashed lines correspond to independent 
emission, quark coalescence, resonance gas and QGP, respectively. The solid circles correspond to random charge assigned to 
each particle and open diamonds are the results from PYTHIA. 



The theoretical prediction 0, for the i'-measure in 
a QGP is = 1. It is 2.9 for a resonance gas 0,01 and 
3.26 in a quark coalescence model [^ITT|. 

Before comparing the data from different experiments 
with the above expectations, one must know how the 
measurements depend on the size of the rapidity window. 
This dependence has been estimated in various models 13, 
H^ . but the results depend strongly on the assumptions 
for the rapidity correlator and the width of acceptance. 
Therefore, a model-independent study of the dependence 
of the fluctuations on the size of the rapidity window in 
over the full rapidity domain is called for. 

In addition, one should test whether the correction fac- 
tors given by Eq. (4) are valid. If correct, a rapidity 
window size scaling should be observed in large rapidity 
windows due to the global charge conservation. More- 
over, the rapidity size for the onset of the scaling will 
offer us a valuable scale for the relaxation time of long- 
range correlation caused by charge conservation 0| . 

Due to the limited acceptance in current heavy ion 
experiments 0, 0, 0| , this kind of study can only be 
performed in hadron-hadron experiments, such as NA22, 
which is equipped with a rapid cycling bubble chamber 
as an active vertex detector and has excellent momentum 
resolution over its full 47r acceptance. 

In this letter, the dependence of the event-by-event 
net charge and charge ratio fluctuations on the size of 
the rapidity window is presented for Tr+p and K+p col- 
lisions at 250 GeV/c. Since no statistically significant 
differences are seen between the results for tt"*" and K+ 
induced reactions, the two data samples are combined for 
the purpose of this analysis. A total of 44 524 non-single- 
diffractive events is obtained after all necessary selections 
as described in detail in ^3 ■ Secondary interactions are 
suppressed by a visual scan and the requirement of charge 
balance, 7 conversion near the vertex by electron identi- 
fication. 

The £'((5)-measure in central rapidity windows \y\ < 
ycut=Sy/2 with Sy = 0.2, 0.5, 0.7, 1.0 to 6 is presented 
in Fig. 1(a), where the open circles are the NA22 data 



TABLE I: D(Q) in different phase-space domains. 



phase-space domain 


RHIC 


NA22 


\y\ < 0.35, Pt > 0.2GeV/c, 
A(j) = 7r/2. 


3.86 ±0.028 
(PHENIX) 


3.896 ± 0.025 


\y\ < 0.5, Pt > O.lGeV/c 


2.8±0.05(cent.) 
3.1 ±0.05(peri.) 
(STAR) 


2.786 ±0.015 



and the open diamonds correspond to the results from 
PYTHIA 5.720 |li| (this convention will be kept in all 
the following figures). In order to compare the results to 
those from STAR and PHENIX, data for the same low 
Pt and azimuthal cuts as STAR 15] and PHENIX [l|| 
are also presented. 

The solid circles correspond to random charge assigned 
to each particle, which indeed gives the value of 4 as ex- 
pected, no matter how small the multiplicity is in very 
narrow rapidity intervals. This shows that the accuracy 
of event-by-event analysis hardly depends on event mul- 
tiplicity and thus can be useful even for low multiplicity 
cases ji^. So, the dependence of the data on centrality 
is not caused by an insufhcient number of particles 20] . 

From Fig. 1(a), it can be seen that with the widening 
of the rapidity window, the NA22 data keep decreasing 
from the value close to 4 (as expected for independent 
emission) to 1 (as expected for a QGP) and even below. 
The loss of small-pt particles and the cut in azimuthal 
angle both enhance the fluctuations. The D{Q) obtained 
under the same cuts as PHENIX [Ti] and STAR ||15| are 
listed in Tab. I. Their values are consistent with ours. 

Taking into account charge conservation in large rapid- 
ity windows and leading-particle stopping, the corrected 
measure D{Q) is presented in Fig. 1(b). The results de- 
crease from about 2.9 (as expected for resonance gas) to 
above 1. The scaling appears when the size of the cen- 
tral rapidity window is larger than \y\ < 2, showing that 
the influence of charge conservation and leading-particle 
stopping have been well eliminated by the factor defined 
in Eq. (4) and the correlation length of charge conser- 
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FIG. 2: The dependence of D(-R+) and D{R-) (upper row) 
and their corrected versions (lower row) on the size of the 
central rapidity window |yj < Sy/2. Open circles are the 
NA22 data and open diamonds are the results from PYTHIA. 
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FIG. 3: The dependence of D{Q), D{R+) and D{R') on 
multiplicity in \y\ < 2 (1st row) and \y\ < 3 (2nd row). Open 
circles are the NA22 data and open diamonds are the results 
from PYTHIA. 



vation is about 4 rapidity units. The corrections reduce 
the measure in small rapidity windows and enhance it in 
large ones. Since the influence of global charge conser- 
vation always enhances the fluctuations, i.e., Cy < 1 in 
Eq. (4) , the suppression in small rapidity windows shows 
that the leading-particle stopping is non-negligible. If 
only the effect of global charge conservation is taken into 
account, as in STAR 15], this will always enhance the 
fluctuations and the scaling behavior disappears. The re- 
sults for such a correction are presented in Fig. 1(c). So, 
the data from both STAR and PHENIX exaggerate the 
fluctuations, the former considering only one correction 
and the latter without corrections at all. 

In Fig. 1(d), D{Q) is presented for different positions 
of a unit-width rapidity window. It is almost a constant 
near that of a resonance gas, showing that the charge 
fluctuations are insensitive to the position of the rapidity 
window and that the local charge is non-equilibrium, as 
pointed out in '2l'|. 

We now turn to a similar study of the charge ratio 
fluctuations. Due to the positive charge of the initial- 
state particles, the average number of positively charged 
particles is higher than that of negatively charged ones. 
Therefore, we present the D-measures in terms of the 
charge ratios i?"*" = n~^/n~ and R~ = n^/n'^ in 
Fig's. 2(a) and (b) separately, where events with = 
and n"*" = have been excluded from the analysis of i?+ 
and respectively. It can be seen from the figures 
that i5(i?+) have much larger values than D{R~). Both 
of them have behavior very different from that of net 
charge fluctuations. 

The charge ratio measures corrected according to 
Eq. (4) is given in Fig's. 2(c) and (d). All points are 
above independent emission and increase rapidly with 
the widening of the central rap idity window, in analogy 
with the model calculation for A- A collisions. These 
results show that the corrections proposed for net charge 



in the observed window as given by Eq. (4) are invalid 
for charge ratio fluctuations. 

It is further interesting to check how these measures do 
in recording the change of charge fluctuations with mul- 
tiplicity in different rapidity windows. This is important, 
in particular, because the even- and odd-multiplicity 
distributions coincide in small rapidity windows, e.g. 
\y\ < 2, while separation of them appears in large win- 
dows, e.g., \y\ < 3. 

The dependence of D{Q), D{R+) and D{R-) on mul- 
tiplicity in two rapidity windows is presented in Fig's. 3. 
The following can be observed: (1) First of all, all plots 
show clear multiplicity dependence, while the results 
from PHENIX [lj| in a small central rapidity window 
show that only D{R) depends on multiplicity, while D{Q) 
is independent of it. (2) For \y\ < 2, the fluctuations of 
even and odd multiplicities in terms of net charge and 
charge ratios coincide within the error bars, consistent 
with the coincidence of even- and odd-multiplicity dis- 
tributions in small rapidity windows. (3) For \y\ < 3, 
the D{Q) separate for even- and odd-multiplicities, con- 
sistent with the separation of even- and odd-multiplicity 
distributions in large rapidity windows. The D{Q) have 
almost equal separation distance for all multiplicities. 
While the D{R'^) and D{R~) are separated differently 
for different multiplicities, with very big errors for odd 
multiplicities, as they could be the combinations of very 
different n+ and n~ . These observations show that D{Q) 
is better in recording the change of charge fluctuations 
with multiplicity in different size of central rapidity win- 
dows. 

The above results can be summarized as follows: 
(1) D{Q), D{R+) and D{R-) depend strongly on the 
size of the central rapidity window. (2) D(Q) eliminates 
the influence of global charge conservation and leading- 
particle stopping. Its scaling behavior is observed when 
the central rapidity window is wider than 4 rapidity units. 
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The same corrections are invalid for charge ratio fluctu- 
ations. (3) D{Q) is insensitive to the position of the 
rapidity bin. (4) D{Q), D[R+) and D{R-) all have clear 
multiplicity dependence. D{Q) has a better record in dis- 
tinguishing the charge fluctuations of even and odd mul- 
tiphcities than D{R+) and D{R-). (5) PYTHIA can re- 
produce almost all the data for charge fluctuations, while 
it fails to describe the transverse momentum fluctuations 
in different central rapidity windows . 

In summary, the dependence of charge fluctuations on 
the size of the rapidity window is presented for the first 
time in the full rapidity domain. The correction factors 
for net charge fluctuations given by Eq. (4) eliminate 
the influence of global charge conservation and leading- 
particle stopping. The latter is non-negligible in small ra- 
pidity windows. Due to the incomplete consideration on 
these two corrections, both STAR and PHENIX exagger- 
ate the fluctuations. The scale of long-range correlations 



caused by charge conservation is about 4 rapidity units 
at y/s — 22 GeV/c. The measure in terms of net charge 
fluctuations is better than that of charge ratio ones. 
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